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Abstract 

The evaluation of loop amplitudes via differential equations and harmonic polylogarithms is discussed at an 
introductory level. The method is based on evolution equations in the masses or in the external kinematical 
invariants and on a proper choice of the basis of the trascendental functions. The presentation is pedagogical 
and goes through specific one- loop and two- loop examples in order to illustrate the general elements and ideas. 
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1 Introduction 



The evaluation of radiative corrections in quantum field theory is notoriously a hard task and various meth- 
ods have been proposed in decades to accomplish it, such as Feynman parameters, dispersion relations, low- 
momentum expansions, etc. In the past few years, a new method has been developed, which is based on (i) the 
reduction of the amplitudes to a minimal set of scalar integrals called master integrals and (ii) their evaluation 
by means of differential equations in the masses or in the external kincmatical invariants; the differential equa- 
tions are then solved using a proper basis of special functions, the harmonic polylogarithms ^ElEIEJE]- Our 
aim is to present a simple introduction to this method. Even though the latter has been used to do multi-loop 
calculations, it can also be used to reproduce standard one-loop results. Indeed, we are going to make use 
mostly of one-loop examples to describe the general elements and ideas. As we shall see, computations which 
are rather laborious with older techniques become much simpler with the differential equations and with the 
harmonic polylogarithms. The plan of this note is the following. 

In sec. |3 we describe the tensor decomposition, i.e. how the evaluation of tensor integrals, coming directly 
from the application of Feynman rules, can be reduced to that of scalar integrals. This step is well known 
to many people and, strictly speaking, it does not belong to the method under discussion; it is included for 
completeness. 

In sec. |3] we describe two widely- used schemes to trasform the dependent scalar amplitudes generated with 
the previous step into a smaller set of linearly-independent ones. We present both one-loop and two-loop 
examples. 

In sec.0]we derive and solve the so-called integration-by-parts identities, which allow to reduce the indepen- 
dent amplitudes to a (much) smaller subset, the so-called master integrals. The two main methods of solutions 
are discussed by means of simple examples. 

In sec. we describe the method of the differential equations to analitically evaluate the master integrals. 
We consider a couple of one-loop examples which exhibit many of the general properties of the method. 

In sec. El we overview the main ideas and results of the harmonic polylogarith theory, including also the 
extension of the basis funcion set to describe amplitudes with threshold at s = 4m 2 . 



2 Feynman diagrams 

The evaluation of virtual corrections to a cross section begins with the application of Feynman rules to the 
relevant diagrams. Delicate points are typically the inclusion of the correct multiplicity factors, the signs of 
fcrmion loops and, whenever gauge interactions are present, a convenient gauge choice. Nowadays this step can 
be done in an automated way [HJ. 

Let us consider for instance the top contribution to Higgs production by gluon fusion, i.e. the process 

9 + 9^ H. (1) 

The Feynman amplitude reads: 

M = e IM (p 1 )e u {p 2 )T^(p 1 ,p 2 ), (2) 

where e^(pi) and e„(p2) are the polarizations of the gluons with momenta p\ = and p\ — and T fJ,v (j>i,p2) is 
the following tensor: 

T^P^-^A^/^Tr 
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where N c is the color factor, g the SU(2) coupling and n the dimension of the space-time. By evalutating the 
trace, the tensor comes out to be of the form: 



d n k 

(2tt)" \{k + pi) 2 - mf] [(fc - p 2 ) 2 - m 2 } [k 2 



,21 ' 



(4) 



where iV 1 " (pi ,p 2 ,k) is a tensor depending on the external momenta p\ and p 2 as well as on the loop momentum 
k: many possibile tensor structures are possible with 3 momenta, 



fT, k»k\ fc"rf, tfk», fc"p£, p»k\ rffi, p^p v 2 



(5) 



The following reduction is convenient. According to relativistic invariance, the tensor can be parametrized as 
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T^(p uP2 ) = rfrfT 1 ( g 2 )+pKT2(g 2 )+^T 3 ( g 2 ) T 4 (g 2 ) + <T(pi •p 2 )T 5 (g 2 



(6) 



where g = px + P2 is the Higgs momentum and t^upa is the antisymmetric tensor. The last form factor Tq is 
related to parity violation of weak interactions. Gauge invariance implies: 



which imply, in turn: 



Pi T M „ = p 2 T MV = 0, 



Ti = T 2 = and T 5 



(7) 



(8) 



The above relations can be used as checks of the computation. The form factors 7j can be derived from the 
tensor by means of projectors. We have for instance: 



T, = 



1 



T 



(n - 2)pi ■ p 2 L" -Pi ' P2 
By applying the projector to both sides of eq. Q and taking it inside the loop integral, one obtains: 

d n k N 5 {k, Pl ,p 2 ) 



T 5 = 



where: 



N 5 {k,pi,p 2 ) = 



(2tt)" [(fc +pi) 2 - m 2 ] [{k - P2 ) 2 - m 2 ] [k 2 - m 2 ] ' 
1 



p 2 Pi 

Pl -P2 



N^(k,p 1 ,p 2 ) 



(n - 2) pi ■ p 2 

We have now a scalar numerator instead of a tensor one, depending only on invariants: 

N 5 (k,pi,p 2 ) = P(k 2 ,k-pi,k-p 2 ), 

where P is a polynomial: 2 



(9) 

(10) 

(11) 
(12) 



P(k 2 ,k -pi,k -p 2 ) = ^2 airs (k 2 ) 1 (k ■ pi) r (k ■ p 2 ) s 

= a 00 o + aioo k 2 + a 01Q k ■ p x + a 0Q1 k ■ p 2 + a 200 (k 2 ) 2 + a uo k 2 k ■ p x + 



(13) 



and ai rs are known constants. n max is the maximum number of invariants and depends on the interaction; 
typically, in one-loop computations, n max = 1, 2. By using eqs. QlUp. I|12|l and (|13|l . we obtain for the form 
factor T5: 



— ^2 a lrs 
l,r.s=0 



d n k 



(k 2 Y(k-pi) r (k- P2 r 



(2tt)™ [(k + pi) 2 - m 2 } [(k - p 2 ) 2 - m 2 } [k 2 - m 



(14) 



2 This is true is all local quantum field theories. 
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3 Independent amplitudes 



As we have seen, the calculations of the radiative corrections to a generic process can be reduced to the evaluation 
of scalar integrals, that in a 2 — > 1 process according to eq. (|14[) are of the form 

d n k (P) 1 (k ■ Pl y (k ■ P2 y 



(2tt)« D 1 D 2 D 3 
where we have defined: 

L>! = k 2 - a, 

D 2 = (k+ Pl ) 2 -a, 

D 3 = (k- P2 ) 2 -a (16) 

with a = m 2 a generic mass squared. 

The above amplitudes are linearly dependent on each other and it is convenient to reduce them to a smaller 
set of linearly independent ones. We express the kinematical invariants in terms of the denominators by means 
of the following formulas — the so-called rotation: 

k 2 = Di+ ft, 
k-Pi = \[D 2 -D 1 ], 

k-P2 = \[-Dz + D x ). (17) 

In general, in one-loop amplitudes, denominators form a basis for the invariants. There are two different schemes 
to implement the rotations given in eqs. 1)17(1 : 

1. auxiliary diagram scheme; 

2. shift scheme. 

1) According to the first method, one uses eqs. I|17fl inside cq. (|14H . to obtain: 

"max 

T 5 = Yl k rs Topo(l-l,l-r,l-s), (18) 

l,r,s=0 

where bi rs are known constants and we have defined: 

/d n k 1 
(2tt)" w or dt ■ (19) 

The amplitudes above constitute the linearly independent ones in this scheme. Let us note that a propagator 
Di originally present in the denominator may be cancelled by a term Di, D 2 , D^, ... in the numerator. In 
diagrammatric language, that means that internal line i is shrunk to a point. The reduction to independent 
amplitudes then generates a pyramid of subdiagrams of the original diagram, in which any subset of internal 
lines is contracted to a point. For clariry's sake, one has to evaluate amplitudes of the form 

1 



Di D 2 D 3 

1 D x D 2 1 D 2 1 D 3 



D 2 D 3 ' D 2 D 3 ' D 2 D 3 ' D.Ds' D,D 3 7 D X D 2 D x £» 2 ' 

J_ Ih. £3 Dl D 2 D 3 J_ Dx D3 J_ 

Di Di' Di' Dx' Di' Di ' "' D 2 ' D 2 D 2 ' ' " D 3 ' V ' 
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2) Let us now consider the alternative scheme for the reduction to independent amplitudes. The shift scheme 
makes use of the substitutions l|17|) only when there is an effective cancellation of denominators. In other words, 
the "allowed" substitutions are: 



Di 
k - pi 
D 1 D 2 
k-p 2 
D\ D 3 



a 

1 1 
~D~x~~D 2 
1 1 

~D~i + D~ 3 



(21) 



This way we are left with independent amplitudes of the following kinds: 



1. amplitudes containing only denominators: 

1 



1 



D 1 D 2 D 3 
1 



DiD 2 



D x D 3 D 2 D 3 
1 1 1 

iV D~2 D~3 



scalar vertex, 
scalar bubbles, 
scalar tadpoles. 



(22) 



The bubble diagrams are obtained shrinking any one of the internal lines to a point, while tadpoles are 
obtained shrinking any pair of internal lines; 



2. bubble diagrams with irreducible numerators: 3 

k-p 2 (k-p 2 ) 2 



DxD 2 



Dx D 2 



k-pi 
Dx D 3 ' 



(k-pxf 
DiD 3 



(23) 



and amplitudes which do not contain anymore Dx, i.e. the denominator containing the loop momentum 
squared k 2 : 



k -px 



k -p 2 



2 \2 



(k 2 ) 



k 2 k ■ px 



k 2 k ■ p 2 



(24) 



D 2 D 3 D 2 D 3 D 2 D 3 D 2 D 3 D 2 D 3 ' D 2 D 3 1 

No simplification is possible for the amplitudes in l|24(l as any cancellation in eq. (|21|l is feasible only if Dx 
is present. We then make a shift of the loop momentum k in order to reproduce a denominator containing 
k 2 , such as for instance: 

k^-k-pi, (25) 

so that 



D 2 
D 3 



Dx, 
Di-- 



(k-px -p 2 f 



(26) 



The shift introduces therefore the new denominator D4, not initially present in the diagram. Since -D4 
contains both px and p 2 , one can express k ■ p 2 in terms of k ■ p\ or vice versa. Let us take the first choice: 



k-p 2 



1 



1 



k-px 



DxD 4 2L>4 2Di DxD 4 
The amplitudes (|24|l are then transformed into amplitudes of the form: 

1 k-px (k-pi) 2 (k-px) 3 



DxD 4 > 



Dx Di 



Dx D 4 



DiD 4 



+ (tadpoles); 



(27) 



(28) 



3 These amplitudes as well as the following ones can be treated with a Passarino-Veltman reduction; we choose to use the general 
method valid also in the multi-loop case. 
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3. Tadpoles, involving: 

• "final" amplitudes of the form: 

k-pi k-p 2 (k-pi) 2 {k-p 2 ) 2 k-p 1 k-p 2 

D 1 ' D x ' D 1 ' £>! ' -Di ' 

• reducible amplitudes of the form: 



jfc 2 


k - pi 


fc -p 2 


(fc 2 ) 2 


k 2 k ■ pi 


k 2 k ■ p 2 




D 2 ' 


D 2 ' 




D 2 ' 


D 2 


k 2 


k-pi 


k-p2 


(fc 2 ) 2 


fc 2 fc • pi 


k 2 k ■ p 2 




D 3 ' 


D 3 ' 


D 3 ' 


£>3 ' 


D 3 



The amplitudes in the first line of the above expression are reduced by means of the shift in (|25() , while 
for the second line we make the shift k — * k + p 2 ; the resulting amplitudes are of the form Ij29(l . 



The same paths of reduction to independent scalar amplitudes can be followed for multiloop corrections. 
As a 2-loop example, let us now consider the light-fermion correction to the process |QJ consisting of a ladder 
diagram. The latter describes a gluon pair converting into a light quark pair which converts in turn into a pair 
of W's or Z's annihilating finally into a Higgs boson. The dependent scalar amplitudes are of the form: 

r d n kx d n k 2 P(k 2 ,k 2 ,ki ■ pi,h ■ p 2 ,k 2 ■ pi,k 2 ■ P2,h ■ k 2 ) , . 

~J (27r)«(27r)« D 1 D 2 D 3 D 4 D 5 D 6 ' 1 ' 

where: 



Dt = kj, 

D 2 = (fci+Pi) 2 , 

D 3 = {k l ~- P 2) 2 , 

Da = k 2 , 

D 5 = {ki + k 2 + pi) 2 - a, 

D 6 = (k 1 + k 2 - P2 ) 2 -a. (32) 



The conversion to independent amplitudes is not straightforward in this case because there are six denominators 
and seven invariants. The denominators then do not form a basis for the invariants, as it happened in the one- 
loop case. 

The solution to this problem, in the auxiliary diagram scheme, is to construct an auxiliary diagram with an 
additional, fictitious denominator linearly independent from the previous ones, such as for instance 

D 7 = (h + k.2) 2 . (33) 

In diagrammatic language, we may say that we have "opened" the Higgs vertex: the auxialiary diagram is a 
planar double box in a forward configuration, i.e. with final momenta equal to the initial ones p\ and p2- After 
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the addition of the auxiliary denominator, the rotation is possible by means of the formulas: 





1.2 
K l 


— ^1 ! 








k 2 


= D 4 , 






h 


■Pi 


= 


Di), 




ki 


■P2 




+ D A ), 




k 2 


■Pi 


- \<*~ 


D 2 + D 5 


-D 7 + a), 


k 2 


■Pi 




D 4 -D 6 


+ D 7 -a), 


ki 


■k 2 


= \(D 7 - 


D x - 





(34) 



The numerator can then be expanded in powers of the denominators as: 

"max 

P(klklk 1 - Pl ,h-p 2 ,k 2 - Pl ,k 2 -p 2 ,k 1 -k 2 )= c hl2hUhl6h D 1 ^ D l » D l i D l l D l i D l i D\\ (35) 

where cz 1 z 2 ; 3 / 4 z 5 ; 6 / 7 are known constants. By inserting the expansion l|35() in eq. (|31|) . one obtains independent 
amplitudes to be computeted, containing formally only denominators, of the form: 

' (36) 



/d n ki d n k 2 
(2tt)" (2tt)" DI 1 



L>™ 2 L»™ 3 L>™ 5 Dq 6 L>™ 7 



with rii < 1 for i = 1 . . . 6 and 717 < 0. Note that the auxialiary denominator D7 appears only in the numerator 
while the standard denominators Di(i = 1 . . . 6) appear both in the denominator and in the numerator. 

As for the case of the auxiliary diagram scheme, also the shift method can be extended to the multi-loop 
case in a straightforward way. 



4 Integration by parts identities 



The virtual corrections to a cross section involve, in general, the evaluation of a large number of independent 
amplitudes — in the case of massive two-loop computations hundreds if not thousands. In the past all these 
amplitudes had to be individually computed It is however possible to reduce by a large amount the number 
of amplitudes to be computed by using integral identities [Hj- In sec. 14.11 we discuss the derivation of the 
identities, while in sees. 14.21 and 14.31 we present two methods for their solution. We work in the auxialiary 
diagram scheme; the discussion in the shift scheme is completely analogous. 



4.1 Derivation of the identities 



Let us begin with the simplest case, that of the one-loop tadpole. According to the divergence theorem: 
where 

L>i = k 2 - a, (38) 



G 



is a sphere of infinite radius in momentum space and ds^ is a surphace element. The flux integral actually 
vanishes only for n\ > n/2, but we will analitically continue eq. I|37[) to all the {n, tlx) space. 

By explicitly performining the derivative and re-expressing the result in terms of independent amplitudes 
by means of the relation (see previous section) 

k 2 = D 1 +a 1 (39) 
we obtain the following integration-by-parts (ibp) identity: 

(n-2n 1 )T(n 1 )-2an 1 T(n 1 + l) = 0, (40) 

where we have defined: 

T(m) - / d n k^. (41) 



By introducing the identity operator I and the plus and minus operators, 
IT(m) = T(m), 

l^K) = T(m±l), (42) 
the ibp identity can be written as: 

[(n-2n 1 )I-2an 1 l + ]T(n 1 ) =0. (43) 



g(n X) n 2 )= / <Tfc nni nW2 , (44) 



Let us now consider as a less trivial case: a bubble with one massive line, 

1 

d';~d. 

where 

£>i = k 2 - a, 

D 2 = (k+p) 2 . (45) 
The integration-by parts identities are derived according to: 

where — fc M or p M . We have therefore two identities for each set of indices (rii, n 2 ): 

(n - m - 2n 2 )I - n x [a (1 + a) + 2"] 1+ = 0, (47) 
(m - n 2 )I + m [a (1 - a?) - 2"] 1+ + n 2 [a (1 + x) + 1"] 2+ = 0, (48) 

where x — —p 2 /a and all the operators are intended to be applied to i?(ni,n 2 ). Three different kinds of 
operators do appear in the identities: 

I, i+, i+j-, (i^j = 1,2). (49) 

The generalitation to multi-loop multi-leg amplitudes is obvious. In the case for instance of the two-loop ladder 
diagram of the previous section, we have: 



- = 0, (50) 



"4 nis n™ 6 n 

^5 u 6 u l 

where j — 1,2 and v = k\, A; 2 , p±, pi is any one of the loop or external momenta. We have eight identities 
for any choice of the indices. Let us note that, in general, the identities are not all independent on each other. 
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4.2 Symbolic solution 



Once the identities have been generated, the next step is to solve them in a convenient way [Hj- Let us begin 
with the tadpole identity 140J1 . We can solve it with respect to the amplitude with the greater index: 

— , / x n — 2ni „ . . 

Tni + 1=- -T(ni) . (51) 

la n\ 

If we assume for instance T(l) to be known, we can determine from the above equation T(2), T(3), T(4), • • • , 
so that we can write: 

T(k) = a(k)T(l), (52) 

where a(k) is a known coefficient and k is an integer. We can also solve eq. ()4(Jf) with respect to the amplitude 
with the smaller index: 

T( ni ) = Ji£_T(n 1 + l). (53) 
n — ln\ 

By setting n\ — we obtain T (0) = 0, and hence T (ni) = for n\ < 0, as well known from elementary 
quantum field theory computations. The conclusion is that the tadpole topology has one master integral, which 
can be taken as T{k) with k a positive integer. 

Let us now solve the bubble identities (|37|) and l|4*5|) . It is convenient to introduce the sum of the indices 

E = m + n 2 . (54) 

The plus operators 1 + and 2 + increase E by one, while the identity I and the plus-minus operators 1+2 - and 
2+1 - keep E unchanged. Let us assume that this topology has one master integral, which we take as 5(1, 1), 
having E = 2 — this will be proved a posteriori. A general amplitude, with ri\ > 1 and n-i > 1, has E > 2. 
That means we have to reduce E by solving the above identities with respect to the plus operators. The first 
equation is solved with respect to 1+: 

l+= n-nr-2n, l _l_ 
ani(l + x) a(l + x) 

With this equation we can shift the first index n\ > 1 down to the value m = 1. Let us remark that it is 
impossible to go further because the coefficients have n\ in the denominator and then become singular. Similarly, 
the second equation can be used to shift the second index ri2 down to one: 

2+ = — 2 - Li+ + -L -1+2- -— rl"2+. (56) 

0,712(1 + X) 77.2(1 + X) 0,712(1 + X) a(l + X) 

Because of the presence of the minus operators, amplitudes with one of the indices equal to zero such as 5(1, 0), 
B(0, 1), B(2, 0), etc., are encountered. These amplitudes have one of the internal lines shrunk to a point and 
are therefore tadpoles, whose reduction has already been discussed. By recursively using eqs. (|55|) and l|56|l we 
can reduce any amplitude B(ni, ri2) with ni > 1 and 712 > 1 to 5(1, 1) + (tadpoles): 

B(m, n 2 ) = c(m, n 2 ) 5(1, 1) + d(m, n 2 ) T(l), (57) 

where c(n 1 ,n 2 ) and d(ni, n 2 ) are known functions. We have thus proved that the one-mass bubble has one 
master integral. 

In some cases, amplitudes of a given topology can be reduced to subtopologies, i.e. to amplitudes with less 
internal lines. Let us consider as specific example a vertex diagram representing the annihilation of two massless 
particles with momenta pi and p 2 into a virtual particle with momentum q = pi + p 2 ■ 

V(n 1 ,n 2 ,n 3 ) = [ d n k— 1 (58) 

J [{k-p 2 y-a\ [(k+piY\ [k 2 \ 
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with pi = and -p\ — 0. There is a massive line "on a side" , connecting one of the massless particles with the 
virtual one. The following two ibp identities are easily derived: 

(n - m -n 2 - 2n 3 )I - m(a + 3")1 + - n 2 2+3~ = 0, (59) 
(-n 2 + n 3 )l + m(-xa - 2" + 3")1 + + n 2 2+3~ - n 3 3+2~ = 0, (60) 

where x = —q 2 /a. By solving the first equation with respect to the 1 + operator and substituting the solution 
into the second equation, one obtains: 

a [{-n 2 + n 3 ) -(n-n 1 -n 2 - 2n 3 )x] I + n x a [-2" + (1 + :r)3~] 1+ + n 2 a(l + x)2+3" - n 3 a3+2~ = 0. (61) 

The above equation does not contain anymore plus operators, bringing unknown amplitudes, but only the 
identity and the plus-minus operators. By setting n\ — n 2 — n 3 — 1 it is immediantely seen that the basic 
amplitude V(l,l, 1) is expressed in terms of amplitudes having one of the indices zero, i.e. of bubbles: we 
succeeded in the above-mentioned reduction. 



4.3 Laporta method 

This method has been originally introduced in |1(J| and has since then been widely used for the evaluation of 
2-loop 3-point and 4-point functions in a variety of mass and kinematical configurations 111! 1121 [5] . The idea 

is that of replacing explicit values for the indices = 1, 0, 1 • • • in the ibp identities. This way a system of 

linear equations is generated, whose unknowns are the amplitude themselves. In the simple case of the tadpole, 
for instance, one generates a system of equations of the form: 



2aT(2) - 


(n 


-2)T(1) 


= o, 


4aT(3) - 


(n 


-4)T(2) 


- o, 


6aT(4) - 


(n 


- 6)T(3) 


- o, 



2kaT(k + l) - (n-2k)T(k) = 0. (62) 

The system is then solved with the method of elimination of variables of Gauss. One has to decide which 
amplitudes have to be solved first. In the above example, one could solve first for T(k), then for T(k — 1), and 
so on. In general, a good criterion is the following |13| : 

• We solve first for the amplitudes with the largest number of denominators. More formally, we define the 
recursive parameter: 

T, 1 =Y / 6{n i ), (63) 

i 

where the step function is defined as 9{u) = 1 if u > and zero otherwise, and we solve first for the 
amplitudes with the greatest Ei; 

• Among the amplitudes with the same number of denominators, i.e. with the same value of Si, we solve 
first for those ones with the greatest sum of the indices of the denominators, 

£ 2 = 5>^K); (64) 

i 

• finally, among the amplitudes with the same values for Ei and E 2 , we solve first for the amplitudes with 
the largest number of D.; in the numerator - in the shift scheme, that is the largest number of irreducible 
numerators: 

£ 3 =^n i 0(-n i ). (65) 
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Many amplitudes have, in general, the same values of (Si, S 2 , S 3 ): the choice of the amplitude can be random. 
Furthermore, a given amplitude appears, in general, in various equations: also the choice of the equation can 
be random. 

According to Gauss method, we proceed with the progressive elimination of variables till all the equations 
have been used. The amplitudes which remain on the r.h.s. at the end are the master integrals. Within the 
scheme (Si, S 2 , S3) presented above, the master integrals typically involve amplitudes with unitary denomina- 
tors (rij = 1) and irreducible numerators (rij = 0,-1, —2 • ■ ■ ). If we exchange S 2 with S3, the master integrals 
typically involve amplitudes with denominators squared. In practise, one usually starts with a small linear sys- 
tem and looks at the master integrals coming from its solution. By enlarging the size of the system, a smaller or 
equal set of master integrals is obtained. The idea of the method is that there is a "critical mass" of equations, 
above which a complete reduction to the master integrals occurs. The reason for this is that, by enlarging the 
system, the number of equations grows faster than the number of unknowns 3 . The main virtue of this method 
is that it can be automated in a rather general way. 



5 The method of differential equations 

Differential equations in the masses or in the external kinematical invariants offer a general method for the 
calculation of master integrals. This method allows in principle to compute any loop amplitude which involves 
more than one scale 4 . In sec. H5.1[) we sketch the derivation of the differential equation for the case of the 
one-mass bubble considered before, while in sees. (|5.2|) and (15.31) we describe the general method to solve the 
equation. 



5.1 Generation of the equation 

Let us begin with perhaps the simplest possible example: the bubble with one massive and one massless line 

1 



where 



B(p 2 ) = J d n k 

d n k = a 2 - n/2 



(k 2 -a){k + vf ' 
d n k 



(66) 



(67) 



«W 2 r(3-n/2)' 

This diagram has a threshold at p 2 = in 2 . To obtain the differential equation, we take a derivative of the master 
integral with respect to the external invariant p 2 using the formula: 5 



d 



dp 2B{p2) 2p* r »dpv 



1 p^B(p 2 ) 



(68) 



The partial derivative is taken inside the integral and produces various scalar amplitudes, which are reduced to 
the master integral itself by means of the methods described in the previous sections. The differential equation 
then closes on the master integral itself: 



1 

x 



1 



1 + x 



B(x; e) + e 



1 + a: 



B(x;e) + (l-e) 



1 



1 + a; 



T(e), 



where 



T(e) = 



(69) 



(70) 



4 Bubble diagrams, vertex diagrams with two external particles on the light-cone, etc., having only massless propagators are then 
excluded. In all these cases, the differential equation gives only a dimensional, trivial information. 
5 We could derive with respect to the mass squared a as well. 



10 



is the tadpole divided by a, 

(71) 

and 

x = - — . (72) 
a 

We have included a minus sign in the definition of x so that the (simpler) euclidean region p 2 < corresponds 
to x > 0. The presence of the threshold in p 2 — a is reflected by the term 1/(1 + x) in the differential equation. 
Let us note that in the above derivation there is nothing specific about the one-loop case so the method extends 
trivially to the multiloop case. 



5.2 Initial conditions 

In order to obtain a unique value for the master integral, an initial condition has to be imposed to the general 
solution of the differential equation. That means we have to know the master integral in a given kinematical 
point x. Let us consider our example. Since B(x; e) is regular for x — * 0, it holds: 

lira x-^B(x; e) = 0. (73) 
x^o ax 

Multiplying both sides of eq. (|69|) by x, taking the limit x — > and using eq. (|73fl . one obtains: 6 

B(x = Q;e) = T(e). (75) 

We have thus obtained the initial condition by studying the master integral close to zero momentum and using 
the differential equaiton itself. 

5.3 Recursive solution in e 

An efficient method to solve the differential equation for the master integral involves the e-expansion of the 
equation itself. Eq. i|()9|) is of the general form: 

e) = >4(x; e) B(x; e) + Q(x; e), (76) 

ax 

where the coefficient of the unknown function is a polynomial of first order in e: 

A(x;e) = A (x) +eA 1 (x), (77) 

with 

Ao(x) = — + —- , 

X 1 + X 

Ai(x) = — — ——- • (78) 
X l + x 

The main point is that A(x; e) does not contain 1/e poles: this is true in general. fl(x;e) is the known term of 
the differential equation and is associated to the tadpole — in general it is related to the subtopologies: 

fi(x;e) = -f2_i(x) +Q (x) + Q.i(x) + Cl 2 (x) + ■■■ , (79) 

£ 



6 In this simple case, the value of the bubble for x = — equivalent to p = — can also be obtained with partial fractioning: 

B(x = 0;e) = / d n k—-^- = - f d"k—^- f d n k\ = - f d n k—^- — = T(e). (74) 

J (k 2 ~ a) k 2 a J k 2 — a a J k 2 a J k 2 — a 

The integral of 1/k 2 vanishes because the integrand is scaleless. 
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where 



0_i(a;) 



flo(x) 

n x (x) 
n 2 {x) = o 



1 1 

X 1 + X 

1 1 



x I + x' 

fii(a!) = 0, 



(80) 



In our example f2(a;; e) has only two non-zero terms; in more complicated cases, il(x; e) contains higher-order 
poles and has an infinite number of positive powers of e. Let us now expand the (unknown) master integral in 
powers of e; since the known term contains at most a simple pole, we expect the same to be true for the MI: 

B(x; e) = -B_x(x) + B Q (x) + eB^x) + e 2 B 2 {x) + ■■■ (81) 
e 

Substituting the expansion (|81|) in eq. H76fl and equating the coefficients of the powers of e, we obtain a series 
of chained differential equations: 

—B- X {x) = A Q {x)B-i{x) + fi_i(x), 
ax 

^Eo(x) - A (x)B (x) + Ax{x)B_ l {x) + Q (x), 

—B x {x) = A (x)B 1 (x) + A 1 (x)B {x) + Q^x), 
ax 

d 



B k (x) = A (x)B k {x) + A l {x)B k _ 1 {x) + (l k (x), 
ax 



(82) 



The first equation, for the coefficient B-i(x) of the simple pole, is the first one to be solved. Once B-\(x) is 
known, we can insert its value in the second equation for Bq(x) and solve for the latter function, and so on. In 
more formal terms, we can redefine the known term as 

£l_i(x) = 

n k (x;B_ 1: --- t B k -i) = A 1 {x)B k _ 1 (x)+n k (x) for k > 0, (83) 
and rewrite the system as: 

d 
i 

d_ 

da 
d 



, B-i(x) = A (x)B- 1 (x) + n_i(ar), 
dx 

. B (x) = A (x)B (x) + Qo(x;B-i), 
dx 



dx 

^B k {x) = A (x)B k (x) + fifcte-B-i,--- 
dx 



(84) 



The system l|84f) is solved with the method of variation of constants of Euler, which we now summarize. Let 
us first consider the associated homogeneous equation for B k (x), i.e. the equation obtained by dropping the 
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known term: 



dx 



u>(x) — Aq{x)oj{x). 



(85) 



An important point is that the above equation is the same for any k. This equation is solved with separation 
of variables: 



ui{x) = exp / A (x) dx 



In our example, 



whose solution is: 



Tx^ 



i 



w(x) = 



1 + x 
1 + x 



w(x), 



x 



where we have taken equal to unity the integration constant. 

The solution of the original, non-homogeneous equation is given by the following integral: 

B k (x)=u(x) ( K(x')n k (x';B^,--- ,B k _{)dx', 



where the kernel K is just the inverse of the homogeneous solution: 

K{x) = -L • 
u>(x) 

The simple pole of the one-mass bubble for instance is given by: 



B_i = c 



l+x 1+x 



dx 



1 



1 



1+x' \ x' 1 + x' 



c-1 



c=l, 



x Jo 

where in the last member we have imposed the initial condition J75}. 

As another example, let us consider the bubble with two equal masses: 

1 



B{x\ e) = / d n k 



[k 2 - a] [(k + pf 



This diagram has a threshold in p 2 = 4m 2 . The differential equation reads: 



1 



1 



2x 2(4 + x) 



B(x- e) 



1 + x 



B(x;e) + (l-e) 



1 



1 



2x 2(4 + x) 



T(e) 



The solution of the associated homogeneous equation in four dimensions, 



1 



1 



2x 2(4 + x) 



<t>{x), 



is: 



<P(x) = 



4 + x 



(86) 

(87) 
(88) 



(89) 



(90) 



(92) 



(93) 



(94) 



(95) 



A first difference with respect to the one- mass case is that the term 1/(1+ x ) is replaced by the term 1/(4 + a;), 
as a consequence of the threshold in p 2 = Am 2 instead of in p 2 = m 2 . Another less trivial difference is that 
semi-integer coefficients appear in the homogeneous differential equation, leading to square roots in the solution 

4>{x). 
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6 Harmonic Polylogarithms 



Let us write explicitly the formal solution of the set of differential equations considered in the previous section: 

f x 

B-i(x) = lu(x) / dx'K(x') Q-i(x'), 

f-X r-X r-X 

B Q (x) = uj(x) / dx'A^x') I dx"K(x")il- 1 (x")+uj(x) / dx'K(x')n Q (x'), 
Bi{x) = lo(x) P dx'Ax{x') I dx"A!(x") f dx"'K(x"') Q-i(x'") 



lu(x) / dx'Ax{x') / dx"K(x")Q (x")+uj(x) / dx'K(x') Qi(a/), 



(96) 

As is clearly seen from the above expressions, the solutions of the differential equations involve repeated integra- 
tions of products of the kernel K{x) and of coefficients of the differential equation itself: Ai(x), f2_i(x), flo(x), 
etc. Natural representations of the solutions seem therefore repeated integrations of K(x) and of the elementary 
functions entering the differential equation under study. The idea behind the Harmonic Polylogarithms (HPLs) 
is simply that of giving a name to such repeated integrations 0]. For the one-mass bubble, for instance, it is 
natural to define: 7 

<?(-!;*) = ' 



l + x' 

g(0;x) = -, 
x 

g(l;x) = -L-. (97) 
1 — x 

The harmonic polylogarithms of weight one are defined as integrals of the above functions: 



H(-l;x) 



r dx' , , 

L TT^ = log(1 + ^ 



f x dx' 

H(0;x) - / —=\og(x), 



H(l;x) = r^ 7 = -log(l-x). (98) 
Jo 1 — x 

Note the slight asymmetry in the lower limit of integration of H(0;x) related to non-integrable singularity of 
l/x in x = 0. Harmonic polylogarithms of higher weight w have the following integral recursive definition: 

H(a,w;x) = [ g{a;x') H(w;x')dx' (99) 
Jo 

for (a, w) (0, W ) and 



H(6 w ;x) = i log w (x). (100) 
wi 

The index a takes the values 0, ±1 and w is a string of w indices, each one taking the values 0, ±1. The vector 
W is a string of w zeroes. 



7 The function 1/(1 — x) is introduced for the closure under the transformation x — + —x. 
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Using the above basis, the bubble with one mass reads: 



B-i = 1, 
B Q = 2- fl-i) // !: 

Bi = 4- (l + -) '2H{-l:.r) + HHl-l:.r) - 2//(-l. -l:.r)] . 



Higher order terms in the e-expansion involve, as expected, harmonic polylogarithms of higher weight. 



(101) 



As far as the bubble with two masses is concerned, the above function set is not sufficient 8 . We add to the 
basis the functions [5]: 



/(-4;z) 
f(-r;x) 



1 



4 + x' 
1 

4-z' 
1 



V^(4+i)' 



/(r;z) - 1 (102) 

^(4 - x) 

The related harmonic polylogarithms of weigth w > 1 are defined analogously to the standard ones. By using 
this extended special function set, the two-mass bubble reads: 

B-i = 1, 



B = 2-4/— — H(-r;x) 



<x 


+ 


4 


4 




+ 


4 


4 



Si - 4- A /— '—[2H(-r;x)-H(-4,-r;x)], 



The lesson is that loop diagrams are represented by complicated, special functions because they involve repeated 
integrations of simple basic functions. In other words, the complexity of the results originates solely from the 
repeated integrations. 
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